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The inspiral and merger of black- hole binary systems are a promising source of gravitational waves 
for the array of advanced interferometric ground-based gravitational- wave detectors currently being 
commissioned. The most effective method to look for a signal with a well understood waveform, such 
as the binary black hole signal, is matched filtering against a library of model waveforms. While 
current model waveforms are comprised solely of the dominant radiation mode, the quadrupole 
mode, it is known that there can be significant power in the higher-order modes for a broad range 
of physically relevant source parameters during the merger of the black holes. The binary black 
hole waveforms produced by numerical relativity are accurate through late inspiral, merger, and 
ringdown and include the higher-order modes. The available numerical-relativity waveforms span 
an increasing portion of the physical parameter space of unequal mass, spin and precession. In 
this paper, we investigate the degree to which gravitational- wave searches could be improved by 
the inclusion of higher modes in the model waveforms, for signals with a variety of initial mass 
ratios and generic spins. Our investigation studies how well the quadrupole-only waveform model 
matches the signal as a function of the inclination and orientation of the source and how the modes 
contribute to the distance reach into the Universe of Advanced LIGO for a fixed set of internal source 
parameters. The mismatch between signals and quadrupole-only waveform can be large, dropping 
below 0.97 for up to 65% of the source-sky for the non-precessing cases we studied, and over a larger 
area in one precessing case. There is a corresponding 30% increase in detection volume that could 
be achieved by adding higher modes to the search; however, this is mitigated by the fact that the 
mismatch is largest for signals which radiate the least energy and to which the search is therefore 
least sensitive. Likewise, the mismatch is largest in directions from the source along which the least 
energy is radiated. 

PACS numbers: 04.25.D-, 04.25.dg, 04.30.Db, 04.80.Nn 



I. INTRODUCTION 

The merger of a binary black hole (BBH) system has 
long been considered a strong source of gravitational 
waves for ground and space based gravitational wave ob- 
servatories. These mergers are characterized by 15 pa- 
rameters, 9 intrinsic to the black-hole systems (2 black- 
hole masses, 2 spin vectors and eccentricity) and 6 ex- 
trinsic to the source (binary orientation vector, sky posi- 
tion and distance). The LIGO and Virgo detectors have 
recently completed a joint run during which inspiral hori- 
zon distances exceeded 40 Mpc [1 and new upper limits 
have been placed on the rates of such events [2 . These 
observatories are currently being upgraded and when the 
new design sensitivities are achieved they will have ranges 
up to ten times greater and hence volumes 1000 times 
greater. By the end of this decade LIGO and Virgo, 
along with GEO, will be joined by KAGRA in Japan 
and possibly the proposed LIGO India, greatly increas- 
ing not only the range of the global network but also the 
ability to recover information about the sources [3j. 

When the theoretical model of the gravitational wave- 
form is well understood, the most effective method to 
search and recover a gravitational wave signal is matched 
filtering against a library of model waveforms called a 
template hank [4 . The ability of such a templated search 
to detect signals is dependent on four factors: 

• The frequency-dependent sensitivity of the de- 



tector. Throughout this paper we use the tar- 
geted aLIGO zero-detuned, high-power [5] sensitiv- 
ity curve. 

• The direction-dependent sensitivity of the detector. 
This is a fixed property of interferometric instru- 
ments and the orientation on the Earth's surface. 
Any one detector will have blind spots, one moti- 
vation for constructing a network of detectors is to 
provide more complete coverage of the sky. We will 
not consider multi-detector searches in this paper. 

• The total energy radiated by the source from the 
time it enters the sensitive band of the detectors. 
This provides an upper limit on the ability to detect 
different signals; a source that radiates less energy 
will be visible out to a smaller distance than one 
that radiates more energy, all other factors being 
equal. 

• The ability of the templates to extract signal power 
from the background noise. 

In this paper we will be concerned with the last two 
points. 

For the BBH systems potentially observable by 
ground-based detectors, astrophysical processes place few 
constraints on the intrinsic physical parameters that 
characterize the emission of radiation from these cata- 
clysmic events, thus placing the burden on source models 
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to cover nearly the full compliment of physical parame- 
ters. Rigorous requirements from matched filtering place 
an additional burden on the source models. In order for 
the model waveforms to match potential signals to within 
a given tolerance, we need not only enough waveforms to 
cover the parameter space but also each waveform must 
represent nature effectually enough to ensure the signal 
does not fall through cracks in the template bank and 
faithfully enough to recover the source parameters. 

One source of mismatch with nature is the truncation 
of the spherical harmonic series in which we have decom- 
posed the model waveform. Current template waveforms 
are only of the dominant, quadrupole mode, although we 
know that generic signals will have many excited har- 
monics present when detected. Fig.Q shows the ratio 
of several non-dominant modes to the dominant mode 
for two non-spinning systems, note that for the system 
where the masses of the component holes are not equal 
the next-to-leading mode is within an order of magnitude 
of the quadrupole mode, suggesting that accounting for 
additional modes may be important for detection, espe- 
cially as the mass-ratio strongly deviates from one and 
generic spins are explored. 

This paper builds on previous work by ourselves and 
other authors. In [6, we conducted a preliminary 
study on higher modes for spinning, equal-mass systems 
comparing numerical relativity waveforms containing the 
largest five harmonics to an equal-mass non-spinning sys- 
tem of just the dominant mode. We found that for low 
spins, the non-spinning dominant mode was an effective 
model waveform. McWilliams et al [8] found that over 
a range of the source orientations, the equal-mass wave- 
form was effective at detecting moderate mass ratios over 
source orientations. Brown et al[9] is exploring the value 
added of higher modes in EOBNR models of unequal- 
mass waveforms. 

In this paper we investigate the degree to which in- 
clusion of additional terms of the spherical harmonic 
series to template waveforms could improve matched- 
filter based searches. We use numerical relativity (NR) 
waveforms as both signal and template, and we consider 
both unequal masses and some generic spins generated 
by the MAYA code. We study how well the quadrupole- 
only model waveform matches the signal as a function 
of the inclination and orientation of the source and de- 
termine how the volume reach of advanced LIGO de- 
pends on the inclusion/exclusion of non-dominant har- 
monics in the model waveforms. We concentrate on sys- 
tem masses greater than IOOMq to give the NR portion 
of the waveform prominence and negating the need for 
post-Newtonian information. Our findings show that for 
non-precessing signals up to 65% of source orientations 
can be missed when using only the quadrupole mode, 
implying a 30% gain in detection volume which could be 
achieved by including higher modes. For our most pre- 
cessing case when using the quadrupole mode only the 
loss of source orientations is 83% and the potential gain 
in volume over which such systems could be detected is 



again 30%. These potential gains in volume are miti- 
gated by the fact that the mismatch is largest for signals 
which radiate the least energy and to which, therefore, 
the search is therefore least sensitive. Likewise, the mis- 
match is largest in directions from the source along which 
the least energy is radiated. Finally, we do a preliminary 
investigation into how the series truncation might impact 
parameter estimation by exploring a potential degener- 
acy between mass and inclination for full waveforms in 
the last section of this paper. 
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FIG. 1: Relative amplitude of higher modes for non-spinning 
Left: q = 1 and Right: q = 4 systems. For the q = 1 system 
the (4,4) and (3,2) modes are about two orders of magnitude 
smaller than the (2,2). All others are less than 10~^. For the 
q = 4 the (3,3) mode is within a factor of 10 of the dominant 
(2,2) mode, and several other modes are within another factor 
of ten. 



We proceed as follows: in 
methodology for matched filtering, and in 
waveforms used in all of our studies. In 



// we introduce our 
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I we con- 
sider various aspects of the overlaps between the dom- 
inant mode and the higher modes. In § [V| we examine 
the volume of the universe accessible to advanced detec- 
tors using quadrupole-only waveforms and hypothetical 
ideal waveforms containing most of the modes, for sev- 
eral cases. We conclude in §|Vl|that the smallest overlaps 



are obtained for systems and source orientations which 
radiate the least total power, and hence have the small- 
est accessible volumes even when an ideal waveform is 
used. In this section we also present a first look at the 
implications of higher modes for parameter estimation. 

Conventions: Throughout this paper we adopt the fol- 
lowing conventions. We denote the Fourier transform of 
a function g{t) with a tilde, as g{f). We characterize 
the mass ratio of a BBH system by q = mi/m2 with 
mi > m2. The relation of the source to the detector 
is specified by five angles. Two (^, ^) place the detector 
in coordinates centered at the source, it is these angles 
in which the decomposition into spherical harmonics is 
performed. Two (^, cp) place the source in the sky of the 
detector. The final angle, tp, determines the relative rota- 
tion between these two coordinate systems, we associate 
iIj with the source because in what follows we will treat 
it similarly to l and (f). We define these angles in fig.([2|. 
The final parameter connecting the source and detector is 
the distance between them, we will be concerned with the 
maximum distance at which the source can be detected 
and will determine this value in what follows. 
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FIG. 2: Definition of angles used in this paper. Left: the 
angles used at the detector, looking at the source. Although 
ip refers to a rotation of the plane containing the source, we 
associate it with the detector because it enters the analysis 
though the antenna pattern. Right: the angles used at the 
source, looking towards the detector. These are the angles in 
which the spherical harmonics are written. 



II. MATCHED-FILTER SEARCHES FOR 
GRAVITATIONAL WAVES 

The response of an interferometric detector is de- 
scribed by an antenna pattern [10 , 
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(1 + cos^ 0) cos 2(p cos 2ijj — cos sin 2(p sin 2?/^ , 
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(1 + cos^ 0) cos 2(p sin 2ijj — cos sin 2(p cos 2?/^ . 



(1) 



Following [TT] we rewrite this in the more convenient form 

F+ = Focos2(V^ + V^o), 

Fx =Fosm2(V^ + V^o), (2) 

where 

Fo = ^((1 + cos2 0) /2)2 cos2 2ip + cos2 sin^ 2cp 

and 



tan2?/^o = 



cos^ 



; tan 2ip . 



(l + cos2 0)/2 

For reference we show the antenna pattern in fig.([3|. 

For gravitational waves, the intrinsic characteristics of 
a source are fully encapsulated in the polarization strains 
h-^ and hx • When an incoming gravitational wave is in- 
cident on the detector the strains give rise to a signal s 
given by 

s{0,ip,L,^,tl;,t) 
= F^{0, ^)h+{L, (/), t) + Fx {0, ^)/ix 0, t) 
= Fo{0,^)h{^,L,^,t), (3) 

where we have used eqn.([2| and defined 

L, (j), t) = COS 2{^l) + ?/^o)^+(^, 0, ^) 

+ sin 2(?/^ + V^o)^x ^, t) 
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FIG. 3: Antenna pattern for an interferometric gravitational- 
wave detector in source-centric coordinates, Lp horizontally 
and vertically. The arms lie along = 7r/2,(^ = 0,7r/2 re- 
spectively. Such a detector is most sensitive to signals directly 
overhead or below, and least sensitive to signals in the plane 
of the arms. The sensitivity drops to zero along the lines be- 
tween the arms, ^ = 7r/2, = zb7r/4 and ^ = 7r/2, = dz37r/4. 



The output of the detector is then 5 + n, where n is 
the noise of the detector. Following standard practice we 
incorporate the noise only as Sn{f) and do not add it to 
the signal in what follows. We will take h in eqn.([3| to 
be the output of a numerical simulation, to be discussed 
in the following section. 

We now briefly review some of the data analysis frame- 
work employed in current LIGO/ Virgo searches, and 
which will be used throughout this paper. An inner prod- 
uct on the space of real, time-dependent waveforms A{t) 
and B{t)^ with respect to a given noise curve described 
by a power spectral density 5'n(/), is 



{A{t)\B{t)) = 4Re df 
Jo 



Mf)B*{f) 

Snif) 



(4) 



In stationary, Gaussian noise, the optimal measure of 
the presence of a gravitational wave signal that matches a 
model waveform, called a template, is the signal-to-noise- 
ratio (SNR) denoted by p, with 



P 



(^+ I ^+) ' 



(5) 



and where we are studying the response of a single de- 
tector to one polarization, typically taken to be We 
note in passing that in a multi-detector search the data 
streams from all instruments will be filtered against the 
same and that the source angles (j) will be the same 
at all detectors. However, the orientations of the differ- 
ent detectors will provide different values of V^, making 
the detectors sensitive to different combinations of the 
polarization. In addition each detector's Fq will have a 



4 



different dependence on 0^(f providing coverage of regions 
of the sky to which any one detector might be insensitive. 

The signal wih arrive at an unknown time which we 
identify as the time of coalescence and denote to. We as- 
sume the template waveform is a good approximation 
to the signal 5, and search for the signal at all times by 
shifting the template. This has the effect in the Fourier 
domain of changing h{f) to exp(— 27rz/to). The sig- 
nal will also have an unknown phase at the time of coa- 
lescence, corresponding to the value of (j) in fig. ([2]), which 
we denote This introduces an additional factor of 
exp(27rz0o)- In practice, this leads the SNR to be evalu- 
ated as 



p{s,h,to) 



I 

Jo 



-27Tifto 



Snif) 



df 



(6) 

where the absolute value removes the dependence on the 
unknown phase. Eqn.([6| may be evaluated by a single 
complex inverse Fourier transform, and the maximization 
over to is then accomplished by finding the maximum 
of the resulting time series. Eqn.(|6| is only an exact 
calculation of the SNR if | ^x) = [E], which is not 
true in general; however, we expect the errors introduced 
by this approximation to be small. 

Note that, by eqn.Q, the dependence on the SNR of 
the detector angles may be factored out in eqn. (|6|. Note 
also that Fo(0,0) = 1. These imply that, given the SNR 
of a signal at 6> = (/:? = 0, we know the SNR of a signal in 
the same orientation at all other sky positions. 

Related to the SNR is the match or overlap obtained 
by normalizing both waveforms 



{s\h. 



max 



{s\h. 



(7) 



The overlap is a measure from to 1 of how well the 
template matches the signal, an overlap of 1 indicates 
that the template is an exact match to the signal and 
anything lower than one is a diminished match. 

Gravitational-wave strain falls off as the the reciprocal 
of the distance between source and detector. It follows 
from eqn.([5| that the SNR falls off in the same way, while 
the normalization removes the distance dependence of 
the template. Henceforth we place the signal s in eqn.([5| 
at 1 Mpc from the detector and denote the resulting SNR 
as piMpc- We also choose a threshold SNR, a value above 
which indicates the presence of a signal in the data. We 
will take this to be 5.5, the threshold used in current 
LIGO/Virgo searches. The choice of this value is moti- 
vated by the behavior of the noise in the detector [13 . 
The distance at which a signal would have an SNR of 5.5 
is then 



PlMpc 

5.5 



(8) 



We now consider two templates, hideai which exactly 
matches the signal and h which in some way approxi- 
mates the signal. We can determine the fraction of the 



available distance that is lost by using the approximate 
template as: 



PiMpc(^)/5.5 



{s\h) 



ideal/ 



{s\h). (9) 



The first equality follows from eqn.(|8]), the second from 

dividing both numerator and denominator by the com- 

1/2 

mon factor {s\s) ' and the third from the fact that when 
the template exactly matches the signal the overlap is 
1. The overlap therefore measures the fraction of the 
SNR lost by using an incorrect template, and equiva- 
lently the fraction of the distance lost. As the universe 
is approximately uniform at distances accessible to even 
initial LIGO [2 , the event rate is approximately equal to 
the cube of the range, although this will also depend on 
the antenna pattern. However, we note that the overlap 
does not give the value of r ideal- As an extreme exam- 
ple, if rideai is sufficiently small that the number of ex- 
pected events per observation time is close to zero, then 
the fractional loss of range implied by a low overlap is 
inconsequential. 



III. THE BINARY BLACK HOLE 
COALESCENCE WAVEFORMS 

This paper uses NR waveforms covering the late in- 
spiral, merger and ringdown for a variety of mass ra- 
tios and spins. All of the NR simulations used in this 
study were produced with GATech's Maya code p^HTO]. 
The Maya code uses the Einstein Toolkit [20 which 
is based on the CACTUS ^ infrastructure and CARPET 
[22] mesh refinement. We use sixth-order spatial finite 
differencing and extract the waveforms at a finite radius 
of 75M, where M is a code unit set to unity and can be 
scaled to any physical mass scale. All grids have 10 levels 
of refinements unless noted below. 

We use 28 simulations in this paper and group them ac- 
cording to their initial parameters in Table [l| Grid details, 
including outer boundary and resolution on the finest are 
also shown. The simulations can be separated into three 
groups: non-spinning, equal-mass with aligned spin, or 
unequal-mass with precessing spin. For the simulations 
with > 4, we used the coordinate-dependent gauge term 
as described in Refs. [22 and [24]. For the g = 10 and 
q — \h simulations, initial parameters in Ref. [25 were 
used. These simulations > 4) have an extra level of 
refinement for 11 levels total, with the exception of g = 6 
and q = 15. These have 10 levels and 12 levels, respect- 
fully. 

The output of all simulations is the Weyl Scalar, 
^4, decomposed into spin- weighted spherical harmonics. 
Simulations are performed in a coordinate system which 
we will denote the source- centric frame, to distinguish 
it from the detector- centric frame we will employ subse- 
quently. See fig.([2| for the definition of the angles used 
in this frame. In terms of these angles the decomposition 
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To quantify this we plot the time series of both SNRs on 
the right-hand side of eqn.(12) in fig.([5|. The two series 
peak at notably different times, and at the peak of the 
h22 series the hother series has dropped by 38% thus we 
can conclude that the non-linearities are important, and 
we cannot use the linear approximation. 



FIG. 4: Overlap of higher modes with 2,2 for Left: q = 1 
and Right: q = 4 systems. In both cases the most significant 
higher modes have poor overlaps with (2,2), suggesting that 
a /i22 will be a poor fit to the full signal in regions dominated 
by these modes. 

is: 

rM^4(^ ^.t)=Y^ -2Yim{i. ^)Cim{t) • (10) 

1,171 

This is related to the strain measured by gravitational- 
wave observatories as 

= ^-2l£m(^,(/>)/iLW. (11) 

£m 

The quadrupole mode is given by (£, |m|) = (2,2). 
Throughout this paper we work in the frequency do- 
main, and therefore avoid the integration to strain since 
h = *4/(-47r2/'). 













Template=/i22 

Template=/i,,/,,,, _ 





























































0.96 0.98 1.00 1.02 1.04 

Time (sec) 



FIG. 5: SNR time series for p(s, /i2,2) and p{s, hothers)- The 
specific behavior will depend on the angles, the values here 
were chosen to illustrate the issue, = 2.36, (p = 2.58,^ = 
1.54, (j) = 5.16. At the time when the /i2,2 series peaks, hothers 
has dropped by 38%. The tension between the modes means 
that the total SNR will be less than the sum of the component 
SNRs. 



IV. OVERLAP 

We start by examining the relative importance of the 
non-dominant modes in a waveform comparison. The full 
waveform involves factors of the spherical harmonics and 
the amplitudes of the modes (see eqn.(lO)). When the 
amplitudes of the higher modes are vanishingly small, 
they can be ignored; however, as we have already noted 
in fig.Q, the relative amplitudes grow in strength with 
mass ratio. 

In fig.Q we plot the overlap of each mode against (2,2) 
individually. If all modes matched well against (2,2) it 
would suggest that a template containing only this mode 
would be a good match to the full signal, regardless of 
the source orientation; however, we find that not to be 
the case. In both the q = I and q = A case, the over- 
lap between (2,2) and the next most dominant modes is 
poor, below 0.6. Furthermore, although the inner prod- 
uct, eqn.Q, and the decomposition into modes, eqn.(lO), 
are themselves linear, the maximization over time and 
phase introduces nonlinearities. In particular, defining 



''Others 



hirai the sum is only linear if the in- 



ner products maximize at the same time. If not, there 
will be a "tension" in the modes and the combined SNR 
will be less than the sum of the individual SNRs, i.e. 



p^{s, h) ^ p^{s, h22) + P^{S, hothers) - 



(12) 



While fig.Q shows that the (2,2) mode is not an ef- 
fective representation of the other modes, how well does 
the (2,2) mode cover the sky of the source? The over- 
lap between the full- mode waveform and the (2,2) mode 
is a function of the angles centered at the source, (^, 0). 
The (2,2)-only template depends on the angles through a 
single factor, 2^22(^7 which is canceled by the normal- 
ization; and, therefore, we simplify the overlap by placing 
this waveform at ^ = ^ = 0. We also place both wave- 
forms optimally in the sky of the detector, at ^ = = 0, 
and choose = 0. We will generalize this momentar- 
ily. Fig.([6| shows the resulting overlaps for five cases, 
the non-spinning q = I and g = 7, and the precessing 
cases from tab.(|l|). At ^ = 0,7r the waveform is domi- 
nated by the (2, 2) modes, the overlap approaches 1.0 at 
these points. Equation ([9| then implies that there is no 
loss of distance incurred by searching with the (2,2)-only 
template for systems that are oriented face-on with re- 
spect to the detector. We can further quantify this by 
determining the faction of surface area over which the 
overlap falls below 0.97%, where this value is motivated 
by the allowed 3% loss of SNR from using a discrete set 
of templates [2 . Tab.([TT|) lists this value for several sim- 
ulations, along with the the average, median and lowest 
overlaps as further measures of the impact of the higher 
modes. 

Figures (5-7) and tab. (O), all tell the same story for 
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409.60 


400 


H09 


15.00 


0.936224 


0.057566 


0.0 




0.0 


(-0.00016052, 0.029072, 


0) 


7.25 


409.60 


800 


SOI 


1.00 


0.453711 


0.453711 


-0.4 




-0.4 


(-0.00079326, 0.092237, 


0) 


11.00 


409.60 


200 


S02 


1.00 


0.453865 


0.453865 


0.4 




0.4 


(-0.00065074, 0.088023, 


0) 


11.00 


409.60 


200 


S03 


1.00 


0.303458 


0.303458 


0.8 




0.8 


(-0.00060332, 0.086010, 


0) 


11.00 


409.60 


200 


POl 


4.00 


0.655334 


0.156900 


(0.6, 0.0, 0.0) 


(-0.6, 0.0, 0.0) 


(0, 0.066502, 0) 




9.00 


409.60 


140 


P02 


4.00 


0.655306 


0.156762 


(0.3, 0.0, 


-0.5) 


(-0.6, 0.0, 0.0) 


(0, 0.068787, 0) 




9.00 


409.60 


140 


P03 


4.00 


0.655306 


0.156764 


(-0.3, 0.0, 


-0.5) 


(-0.6, 0.0, 0.0) 


(0, 0.068758, 0) 




9.00 


409.60 


140 



TABLE I: Simulations Used: The 28 simulations' initial parameters and grid structures are listed. The table is split into 
three groups: non-spinning, equal-mass with spin, and precessing spins. The table contains q = m+/m-, the bare puncture 
masses mh+/M and mfe_/M, the non-dimensional spins, Xi — Si/'^l^ the initial momentum, p+/M, the initial separation, 
d/M, the outer boundary, Rb/M, and the resolution on the finest refinement level M/hfine- If only one spin value is listed, 
the spin is aligned with the initial angular momentum. 



a single detector when the intrinsic parameters are kept 
fixed to the signal: the q=l case is well served with a 
(2,2)-only waveform over all source angles. The higher 
the mass ratio, the worse a (2,2)-only waveform does 
in matching the signal, and this fraction of angles over 
which the match does poorly increases. Furthermore, a 
precessing system is badly served by a (2,2) waveform. 
We will explore this matter further in future work. 

We now generalize the previous results to include other 
values of the detector-centric angles (^, ip) and ?/^. Con- 
sider two templates: /122 which, as in current searches, 
contains only the (2, 2) mode of the NR waveform opti- 
mally oriented (^ = (/9 = ^ = (/) = ?/^ = 0), and a per- 
fect template hideai which exactly matches the signal. In 
we show the overlap between the signal and /122 for 
several non-spinning systems. Each colored line on the 
graph represents a system mass ratio, moving along the 
line gives different system masses. As we move from top 



to bottom, we are moving from q=l to q=15. The differ- 
ence in colors along the line give the overlap value. The 
plot shows that for higher mass ratios the total power is 
distributed into the higher modes and the match drops 
accordingly. This is consistent with |8l l9j . 

Now consider the q = A non-spinning system, scaled to 
to IOOM0 and placed at a distance of IGpc from the de- 
tector, and examine the overlap between the signal and 
both templates. We randomly choose values for all an- 
gles and plot results with respect to ^, which has the most 
significant dependence. The results are shown in fig.Q, 
which illustrates that at ^ = 0, tt the variation of the addi- 
tional angles do not affect the overlap, while the spread in 
results widens towards i = 7r/2. This again shows that 
the (2,2) mode only captures a face-on source orienta- 
tion and misses the source as its inclination increases to- 
ward the edge-on case. This would imply that the higher 
modes are essential for detecting non-optimally oriented 



7 







0.96 
0.94 
0.92| 
0.9oJ 



0.98 
0.96 
0.94 Q 
0.92 S 



ID 


(I 


a 


% of area 


Average 


Median 


Minimum 








> 0.97 








HOI 


1 





100 


0.997 


0.998 


0.995 


H03 


3 





51 


0.955 


0.951 


0.918 


H04 


4 





43 


0.937 


0.931 


0.885 


H05 


5 





40 


0.927 


0.920 


0.868 


H06 


6 





37 


0.916 


0.907 


0.847 


HOT 


7 





36 


0.907 


0.898 


0.840 


H08 


10 





36 


0.903 


0.892 


0.826 


H09 


15 





35 


0.897 


0.886 


0.817 


SOI 


1 


-0.4 


100 


0.997 


0.997 


0.993 


S02 


1 


0.4 


100 


0.997 


0.997 


0.994 


S03 


1 


0.8 


100 


0.997 


0.997 


0.994 


POl 


4 


0.6 (90°) 


13 


0.883 


0.889 


0.741 


P02 


4 


0.6 (150°) 


41 


0.938 


0.939 


0.852 


P03 


4 


0.6 (210°) 


28 


0.933 


0.942 


0.816 



TABLE H: Summary values of the overlaps between the (2,2) 
mode and the full template as a function of the orientation 
angles (^, 0). Names in parenthesis refer to tab.(|l]). Note that 
the POl precessing system has lower overlaps, and a smaller 
fraction of overlaps greater than 0.97, then the other systems. 



FIG. 6: Overlaps in source-centric coordinates, <j) horizontally 
and L vertically, between the complete waveform and the (2, 2) 
mode for Top: the non-spinning q — 1 and ^ = 4, Middle 
the precessing POl and P02 and Bottom: the precessing P03 
signals from tab.(|l]). The general features of the non-spinning 
images are representative of all mass ratios and (anti-) aligned 
spin systems; overlaps are 1.0 at ^ = 0, tt where the full signal 
reduces to the (2,2) mode, and are lowest at ^ = 7r/2. There 
is more interesting structure in the precessing cases. 



signals, but how far away can a single detector see these 
cases? We quantify how important the modes will be in 
terms of SNR and volume reach in the next section. 



SNR AND VOLUME 



^ 60 

e" 




50 100 150 200 250 



As noted at the end of § [77j the overlap is equal to 
the fractional loss in distance to which a signal can be 
detected, but this value should be viewed in light of the 
maximum possible distance. This maximum distance de- 
pends on three factors: (1) the total energy radiated by 
the source, (2) the ability of the template to extract en- 
ergy of the signal from the background noise and (3) the 
location of the source in the sky of the detector. For 
example, in the plane of the detector along the lines 45 
degrees to the arms, the response goes to zero. Along 
these lines the loss in range implied by a low overlap is 
irrelevant for a single detector. In this section we con- 
sider the accessible distances, noting the influence of all 
three factors. 

We start with fig.(|9|, which shows the radiated energy 
and distances accessible using the hideai templates, as 



FIG. 7: Overlaps between the complete waveform and the 
(2,2) mode for non-spinning waveforms with mass ratios from 
1 to 15, with all angles and total mass chosen randomly. At 
higher mass ratio more of the total power is distributed into 
the higher modes and the match drops accordingly. 



a function of the source orientation. As expected, the 
range tends to be lowest where the least power is ra- 
diated, although the energy and distance plots are not 
identical due to weighting by the noise curve. The en- 
ergy, and hence distance, plots have the same general 
shape as those corresponding in fig.([6|, indicating that 
the overlaps between the signal and /122 are lowest at 
orientations where the energy and distance reach of the 
ideal template are also lowest. This is due to the fact that 



FIG. 8: Left: The overlaps obtained using both templates. 
Since hideai = s the overlap is 1. Right: The ratio of the 
overlaps. This is identical to the ratio of SNRs as the addi- 
tional factors of (s | s) cancel. 



the higher modes not only have poor matches with (2, 2), 
as shown in fig.Q, but they also contain less power, as 
shown in fig.Q. Fig.(|9]) shows that orientations where 
the higher modes dominate have both low matches with 
h22 and lower ranges. This indicates that the fractional 
loss in distance incurred by using the incorrect template 
is greatest where the best possible range is smallest. 

Finally, in order to characterize the performance of 
different templates by a single number with physical sig- 
nificance we calculate the spatial volume to which the 
search is sensitive. The distance to which a signal can 
be seen depends on all five angles, but from eqn.([3| and 



the comments at the end of § // the dependence on the 



detector-centric angles may be factored out 

R{0, (/p, L, (/), ^) = Fo{0, ^)R{l, (/), ^) (13) 

Since there is no preferred orientation we define an 
average visibility range^ by averaging the distances 
over the orientation angles i^(p^^lj: 



1 ^ 



p{s{ii,(l)i,il)i),h) 
5.5 



(14) 



We evaluate this average by choosing random values 
for ^, cos(0), uniform in (0, 27r), (— 1, 1), (0, 27r) respec- 
tively. 

The average visibility distance as a function of the 
detector-centric angles is therefore 



Ri0,ip)=RFoi0,ip) 



(15) 



and the volume of the Universe to which a given template 
is sensitive is therefore 



V 



dip sm{0)dO / r 
Jo Jo Jo 



I p2,7T pTT 

- dip sm{0)deR^{0,ip) 
3 Jo Jo 

dip sm{0)dOF^{0,ip). (16) 
Jo Jo 



3 
3 



The remaining integral may be done numerically, yielding 
a value ~ 3.687. 




1.20 
1.05 J 



0.75 s 



0.60^" 
0.45 2 

0.30 










FIG. 9: Radiated energy and distances to which signals are 
visible using the optimal template, in source-centric coordi- 
nates (j) horizontally and t vertically. Top to bottom: q — 1, 
Q = 7, and the precessing POl, P02 and P03 systems. Note 
the structure is similar to the overlaps between the full signal 
and the /i22 template, fig.Jol. 



The volumes for different waveforms, using the /122 and 



hideai templates are summarized in table IIT The trend 
is for lower mass ratios and higher aligned spins to cor- 
respond to both larger absolute volumes and smaller rel- 
ative differences by including higher modes in the tem- 
plate. The larger volumes correspond directly to the in- 
creased total energy radiated by such systems, which is 
shown in fig. (10). 



Finally, as another way of quantifying the difference 
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10^ 



10^ 



10^ 




).00 0.01 0.02 
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0.07 0.08 0.09 




Distance (Gpc) 



FIG. 10: Correlation between the total energy radiated from 
r = 75M by the systems in table |III| and the accessible vol- 
umes using the /i22 and hideai templates. Circles are non- 
spinning systems, squares are spinning but non-precessing 
systems, and triangles are precessing systems. The P02 and 
P03 systems have close to identical values of Erad and vol- 
umes, these points therefore lie on top of each other. The /i22 
template gives a notably smaller fraction of the volume for 
the POl system than for any other, this corresponds directly 
to the lower overlap noted in tab.([TT|. 



between the templates, in fig. (11) we show histograms of 



the visibility ranges over the complete set of orientations 
dit 6 = (p = 0. Using hideai shifts the ranges from lower to 
higher values somewhat, but does not increase the max- 
imum distance, which occurs for face-on systems which 
are dominated by (2,2). 

These results include three precessing g = 4, a = 0.6 
systems. In all cases the accessible volume is less than 
that for the q = 4 non-spinning system. As might be 
expected from the non-precessing cases the volume de- 
creases as the spin becomes anti-aligned with the angular 
momentum and less total energy is radiated. However, 
at least for the systems considered here, this dependence 
becomes smaller than our uncertainties when the angle 
between the orbital angular momentum and the spin of 
the larger hold exceeds 150 degrees. 



FIG. 11: Histograms showing the distributions of distances 
using both templates for the q = 4 system. Using hideai shifts 
points from lower distances to higher, but does not increase 
the maximum range. 






FIG. 12: Histograms showing variation in distance along 
the = cp = sky direction for Left q=l and Right q=7 
systems. For q = I the mean is 30.80, corresponding to a 
distance of 5.6 Gpc, and the standard deviation is 0.17. For 
q = 7 the mean is 13.79, corresponding to 2.51 Gpc, and the 
standard deviation is 0.08. 



3Sr/r. The error for the results in tab. (Ill) is then on 



order 1.5%. There are also uncertainties associated with 
the choice of extraction radius and resolution. We show 
the volumes obtained using the q = 4 systems and hideai 
template for several value of both parameters in tab. (IV). 
The variation is on the order of 1.5%, and our two sources 
of uncertainty are comparable, and small enough that 
they do not effect our conclusions. 



A. Error analysis 



VI. CONCLUSIONS 



Because we choose random values in evaluating the 
average eqn.( 14) we are able to determine the error in the 



results as the standard deviation between several runs. 
Due to the computational expense of complete runs we 
instead estimate this by choosing one sky position. We 
show the SNR histograms obtained by 900 runs of ^ = 
(f = 7t/3 for two waveforms in fig. (12). In both cases the 



error is on order of 0.5%. Since V 



"3 



and r has an error 



Jr, then SV = y {{dV/dr)Sr) . Here we have 5V/V = 



As can be seen from table HI there are two conflicting 
trends as the mass ratio increases. As the total radi- 
ated energy is reduced, the volume drops. Conversely, 
as the fraction of this energy is distributed into higher 
modes the benefit gained by using the ideal template 
increases. The energy radiated, and hence volume, in- 
crease with spin. Together, these results imply a strong 
bias towards the detection of equal-mass, aligned- spin 
systems when averaged over the sky. This conclusion is 
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ID 


q. 


a 


Volume using 


Ravg 


using 


Volume using 


Ravg using 








h22 (Gpc^) 


h22 


(Gpc) 


hideal (GpC^) 


hideal (Gpc) 


HOI 


1 


0.0 


217 




3.3 


218 


3.4 


H03 


3 


0.0 


91 




2.5 


102 


2.6 


H04 


4 


0.0 


57 




2.2 


68 


2.3 


H05 


5 


0.0 


39 




1.9 


47 


2.0 


HOG 


6 


0.0 


27 




1.7 


34 


1.8 


HOT 


7 


0.0 


19 




1.5 


25 


1.6 


H08 


10 


0.0 


9.3 




1.2 


12 


1.3 


H09 


15 


0.0 


3.3 




0.8 


4.3 


0.9 


SOI 


1 


-0.4 


165 




3.1 


166 


3.1 


S02 


1 


0.4 


313 




3.8 


315 


3.8 


803 


1 


0.8 


458 




4.3 


461 


4.3 


POl 


4 


0.6 ( 90°) 


41 




1.9 


55 


2.1 


P02 


4 


0.6 (150°) 


33 




1.7 


39 


1.9 


P03 


4 


0.6 (210°) 


33 




1.8 


39 


1.9 



TABLE HI: Sensitivity volumes and average distances achievable using both templates. ID values correspond to tab.Q Angles 
following spin magnitude indicate the initial angle of the spin vector of the larger hole in the x, z plane, such systems exhibit 
precession. Spins not followed by an angle indicate the spins are (anti) aligned with the orbital angular momentum and the 
system does not precess. Volumes are reduced with increased q and anti-aligned spins, and increased with align spins due to 
total power radiated in-band. For higher q the use of the ideal template expands the volume by up to 30% for the systems 
considered here, although the fractional improvement is greatest for the systems where the volume accessible with hideal is 
smallest. 



Extraction r. 


Resolution 


M/160 


M/180 


M/200 


60 M 




69.59 


69.62 


69.64 


75 M 




69.12 


69.14 


69.16 


100 M 




68.57 


68.59 


68.61 



TABLE IV: Volumes obtained using the q — 4 system and 
hideal template for various extraction radii and simulation 
resolutions. All values are in Gpc^. All of these runs used the 
same set of points. There is a general trend downward with 
decreased resolution and increased extraction radius. The 
latter effect is due to the fact that the late inspiral, merger 
and ringdown portions of the waveform get smaller as r ^ oo. 
Although the inspiral portion actually increases as r ^ oo, 
since the majority of the power radiated is in the last orbits 
and merger the volume decreases. As the variation is small 
we expect the difference from the true value to be small as 
well. 



consistent with [26l [27] , while adding the fact that the 
inclusion of higher modes is not important for detecting 
these systems. We expect that a search using (2,2) IMR- 
PhenB aligned-spin templates will perform well, this will 
be tested as part of the ongoing NINJA2 project [2F. 

For non-spinning systems with ^ ^ 3 and the mildly 
precessing systems considered here, the inclusion of 
higher modes in the template can improve the volume 
reach of the single detector. Whether or not this trans- 
lates to an increase in detection rate depends on the un- 
known underlying rates of such systems. Put another 



way, the inclusion of higher modes in templates will al- 
low the advanced detector network to better measure or 
bound these unknown rates. 

There are, however, some caveats. First, we stress that 
the template used for the rightmost column of table III 
exactly matches the signal, that is, it assumes we ex- 
actly know the signal for which we are looking in ad- 
vance. To the extent that matched filtering is the optimal 
detection statistic any approximate inclusion of higher 
mode information will of necessity do worse. Further- 
more, there are potential downsides to including higher 
modes in the templates. Such an addition would re- 
quire increasing the number of templates. This entails a 
corresponding increase in the computational cost of the 
search. In addition, these additional templates may re- 
spond to glitches in the detector, raising the number of 
"background" events and increasing the SNR at which a 
signal would need to be observed in order to confidently 
claim a detection. Concerns such as this lead to changing 
the mass range in the S6 search from 35Mq to 25M0 - 
the templates at the higher mass end produced sufficient 
numbers of background triggers to impair the ability to 
detect lower-mass systems [1 . It would be undesirable 
to allow a search for systems to which the detector net- 
work is comparatively insensitive to impact the ability 
to detect equal-mass and aligned-spin systems. We also 
note that, at present, it is not known how to construct a 
template bank of precessing signals. Further studies are 
needed to determine the right strategy for detecting both 
mildly and heavily precessing systems. 

We have not yet considered spinning systems with 
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g > 1. Such simulations are available for spins up to 0.6 
and mass ratios up to 7, however, we defer their anal- 
ysis to future work. For spins aligned with the angu- 
lar momentum the volumes accessible will certainly be 
larger than the non-spinning counterparts. It is possible 
that the dependence on higher modes will be preserved in 
these cases, leading to a potentially large volume increase 
by using templates that include higher modes. 

We have so far considered only a single detector. Ad- 
ditional detectors will provide better sky coverage, ef- 
fectively increasing the value of the integral in eqn.(16). 



Furthermore, as noted in the introduction, detectors ori- 
ented differently are sensitive to different polarizations, it 
is therefore conceivable that the inclusion of higher modes 
in templates would have more impact on the range of the 
network as a whole than on any one detector. We have 
also not considered other aspects of the full search, such 
as signal-based vetoes. The effect of such vetoes is being 
studied in [29j. 

One important aspect of gravitational- wave detection 
we have also not considered is the fact that the data is 
filtered against a bank of templates with different pa- 
rameters. For the initial detection it is acceptable for 
the signal to be picked up by a template with the wrong 
parameters; once the detection has been confirmed more 
computationally expensive parameter estimation codes 
can be run. While this freedom can not raise the volume 
accessible to hideaU as it is already a perfect match to 
the signal, it is quite possible that maximization over a 
bank of /122 templates will lead to larger average SNRs 
and hence volumes. In this case the fractional gain by 
going to an approximation of hideai may be even smaller. 

This last point leads to the question of the importance 
of higher modes in parameter estimation. We expect 
higher modes to be important here; as a simple example 
the difference between a signal dX i = = ii/A and one at 
i = il) = {) is entirely encapsulated in the mode content. 
We expect that there are degeneracies between the orien- 
tation parameters and intrinsic parameters, we intend to 
investigate this further in subsequent studies. However 
we present a preliminary result in fig. (13), which shows 



that {s{i^ij)) \ h22) can be increased by maximizing over 
the mass M of the template, at the cost of misestimating 
the mass. The increase in overlap is most pronounced 
at i = 7r/2, where the higher modes are most significant. 
Correspondingly the mass which maximizes the overlap 
deviates the most from the true value at this point. This 
suggests a degeneracy between mass and higher mode 
content. One possible explanation is that the higher 
modes contain more power at higher frequencies, as do 
lower-mass systems. We will explore this possibility in 
our follow-up studies. 
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